/J 

'J 


f 


00 

C3 

t- 


£ 

s 

E-i 

< 

O 

g 


oo 

ro 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


A 


D ! 

Si 

-C 1 

oj 

tn 3 


TECHNICAL NOTE 


No. 1723 


A THEORETICAL INVESTIGATION OF THE EFFECT OF YAWING MOMENT 

DUE TO ROLLING ON LATERAL OSCILLATORY STABILITY 

By Joseph L. Johnson and Leonard Sternfield 

Langley Aeronautical Laboratory 
Langley Field, Va. 




Washington 
October 1948 


7" r— <** n P r 


d i j 


> / -0 


TECH LIBRARY KAFB, NM 



NATIONAL ADVISOR! COMMITTEE FOB AERONAUTICS 


QS 


IMM'l 


TECHNICAL NOTE NO. 1723 


A THEORETICAL INVESTIGATION OF THE EFFECT OF YAWING MOMENT 
DOE TO R OLLIN G ON LATERAL OSCILLATORY STABILITY 
By Joseph. L. Johnson and Leonard Sternfield 

SUMMARY 


The effect an lateral oscillatory stability of varying the value of 
the stability derivative C_ ( yawl ng-mament coefficient due to rolling— 

“p 

angular—' velocity factor) over a wide range of negative and positive values 
has been determined theoretically for a sweptback fighter airplane with 
loading conditions simulating an airplane with wing-tip fuel tanks and 
without wing-tip fuel tanks. Increasing C^ in a positive direction 

was shown to increase the damping of the short-period oscillation. For 
any given value of the addition of the wing-tip t ank s, which 

involved shifting fuel from the fuselage to the wing tip, decreased the 
damping of the short-period oscillation. 


INTRODUCTION 


Recent theoretical and experimental investigations (reference l) 
have Indicated that the stability derivative (yawing-moment coeffi- 

cient due to rolling-angular-velocity factor) may be greatly affected by 
sweep, aspect ratio, taper ratio, and, lift coefficient. In fact, the 
value of Cj^, which is normally negative in sign and proportional to 

lift coefficient for straight wings, has been found to vary from -0.10 
to 0.30 over the lift range for highly swept wings. Calculations of 
reference 2 showed that even rnnal 1 variations in the value of C n 

(—0.02 to 0.02) might have an appreciable effect on the neutral— 
oscillatory— stability boundary. In order to determine the effect on 
the oscillatory stability of extending the range of the variation of C^, 

stability calculations were made in which C^ was varied systematically 
over a wide range of positive and negative values (—0.10 to 0.40). 


The calculations were made for a sweptback, jet-propelled fighter- 
type airplane for two loading conditions. The loading conditions assumed 
were for the airplane with wing-tip fuel tanks an d without wing-tip 
tanks. For the condition without wing-tip tanks an appr oxi m a tely equal 
quantity of fu^l was assumed to be carried near the center of gravity . 
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SYMBOLS AHD COEFFIdEHTS 


•wing area, square feet 
mean aerodynamic chord, feet 
airspeed, feet per second 
•wing span, feet 

dynamic pressure, pounds par square foot 
air density, slugs per cuMc foot 
■weight, pounds 

acceleration of gravity, feet per second per second 
mass, slugs (W/g) 

relative density factor "based an wing span (m/pSb) 

angle of attack of reference axis (fig. l), degrees 

angle of attack of principal longitudinal axis of airplane, 
positive when principal axis is above flight path at the 
nose (fig. l), degrees 

angle between reference axis and principal axis, positive when 
reference axis is above principal axis at the noBe (fig. l), 
degrees 

angle between reference a xi s and horizontal axis, positive 
when reference axis is above horizontal axis at the nose 
(fig. l), degrees 

angle of flight to horizontal axis, positive in a climb 
(fig. l), degrees 

angle of yaw, degrees or radians 

angle of sideslip, degrees or radians 

angle of bank, radians 

Kouth*s discriminant (B = BCD - AD 2 - b2e where A, B, C, 

D, and E are constants representing coefficients of the 
lateral stability equation) 

radius of gyration about principal longitudinal axis, feet 
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radius of gyration about principal vertical axis, feet 


nondimensional radius of gyration about principal longitudinal 
axis 

-n nnrH-nipmH ln na.1 radius of gyration about principal vertical 
axis j b^ 

•n onfl-i-nipinpi ^ nnal radius of gyration about longitudinal stability 
axis 2 cos 2 t] + E^^sin 2 ^ 

n nrifilTnAnfl l anal radius of gyration about vertical stability 
axis ^j/Kz o 2 cos 2 t| + Ky^sin 2 !}^ 

nondimensional product-of— inertia parameter 
(L, 2 ~ cos t] sin 

lift coefficient (lift/qS) 

yawing-moment coefficient (Tawing moment / qSb ) 

-ro l 1 Ing -moment coefficient (Boiling moment /qSb ) 

lateral— force coefficient (lateral force /qS) 

rate of change of lateral— force coefficient with angle of 
sideslip, per degree or per radian, as specified (SCy/dp) 

rate of change of yaw in g-moment coefficient with angle of 

sideslip, per degree or per radian, as specified (BC^/dp) 

rate of ch ang e of rolling-moment coefficient with angle of 
sideslip, per degree or per radian, as specified (dC^/dp) 

rate of c hang e of lateral— force coefficient with rolling- 

hoj \ 


angular— velocity factor, per radian 


a ^ 

° 27 


rate of change of rolling-moment coefficient with rolling- 

faj \ 

angular— velocity factor, per radian / \ 
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rate of change of yawing-moment coefficient with rolling— 
angular-velocity factor, per radian 




rate of change of rolling-moment coefficient with yawing— 


angular— veloc ity f ant or, per radian 


rate of change of yawin g-moment coefficient with yawing— 

'a9a‘ 

angular— Telocity factor, per radian [ — — 

27 / 

rate of change of lateral— force coefficient with yawing— 
angular-velocity factor, per radian 



tail length (distance from, center of gravity to rudder hinge 
line), feet 

height of center of pressure of vertical tail ah ore fuselage 
axis, feet 

rolling-angular velocity, radians per second 
yawing-angular velocity, radians per second 
differential operator (d/ds^) 
distance along flight path, spans (Yt/b) 
complex root of stability equation (c ± id) 
time, seconds 

period of osc illat ion, seconds 

t ime for amplitude of oscillation to change by factor of 2 
(positive value indicates a decrease to half— amplitude, 
negative value indicates an Inc rease to double amplitude) 
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EQUATIONS OF MOTION 


The nondimensional lateral equations of motion (reference 2), 
referred to a stability-axes system (fig. 2), are: 

In. roll — 


2 ^x 2B b 2 0 + + I + I 


In yaw — 

* Sxsft)^) - V + i + I VW 


In sideslip — 

2 M^( I> bP + + \ ^Tp-^b^ + + 2 ^j^b^ + ( C L ^ an 

Asij Asb Asb 

When (p Q e is substituted for <p, ^ Q e for y, and P D e 

for P in the equations written in determinant form, X must be a 
root of the stability equation 

AA^ + BA 3 + CA 2 + DX + E = 0 (l) 


where 

A = 8pb 3 (% 2 K Z 2 ~ %Z 2 ) 

B = -2Hb 2 (2Kx^z 2 % + + Ez 2 % - 2KSZ 2 % ~ ^XZ^ - %*;%) 

C = Mb(% 2 Qn r °yp + %ib% 2c np + E Z 2c I p °rp + 2 C n r Cz p ~ K XZ C I r °rg 
~ ^^XZ^p “ ^nX-TZpYp ~ \ C n p G Z r + ^Z^ 0 ^ “ p c Zp 


K X 2c T r C np + %Z<Vi p ) 
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d -Kvp 

~ 2 ^tPlpL G l^ 

+ 1 %%% 
E - I °L(Ca/lp - 


- ^ C I ? c np + ^ C B 5 G l I . 0 Yp + ^ C n p C Ip + 2 ^®L%Z°np 

- tan 7 + Zv&jzpi Cfc tan 7 

- i %%%,. - 1 C l T %% * fW, 

°vS) + 1 °l tan 7 (%% - %%) 


The damping and. period of the lateral oscillation in seconds are 
given respectively by the equations ^ 1/2 = ~ ""TT^y an( ^ P = ~ ^ 


•where c and d are the real and imaginary parts of the complex 
root of the stability equation. 


The necessary and, sufficient conditions for neutral oscillatory 
stability as shown in reference 3 are that the coefficients of the 
stability equation satisfy Bouthis discriminant set equal to zero 


E = BCD - AD 2 - B 2 E = 0 


and that the coefficients 33 and 1 have the same sign. In general, 
the sign of the B coefficient is determined by the factors -Cy^, —C^, 

-C7_ which appear in the predominant terms of B. Thus, B is 


and 


'ip 


positive in the usual case in which there is positive -weathercock stability 
and positive damping in roll. Hence the coefficient D must be positive 
if E = 0 Is a neutral— o s c illat ory-st ability boundary. 


CALCUIATIONS 


Calculations were made of the neutral— oscillatory— stability boundary 
for the two different loading conditions representing the airplane with 
wing-tip fuel tanks and, the airplane without wing-tip fuel tanks. These 
two loading conditions simulate a change in the fuel load from, a belly 
fuel t ank to wing-tip tanks and, in effect, vary the radii of gyration 
in roll and yaw while the airplane weight is maintained approximately 
the same. 

The roots of the stability equation were also computed for several 
combinations of Cq , Cjg, and C^ to determine the period and time 
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to damp to one-half amplitude of the oscillatory and aperiodic modes of 
motion. 

The aerodynamic and masB characteristics assumed for use in the 

calculations are presented in table I. Yalues of Cy,., 

^P(tail off) 

and Cy -were estimated from force tests made on the 

P(tail off) 

Langley free— flight-tunnel six-component balance (reference 4) on a 
model representing the type of airplane assumed in the calculations. 

The tail-off values of CL.Ct.Cv.C, , and. Cy were estimated 

•“r V x r V X P 

from data taken In the Langley stability tunnel. The contributions of 
the tail to the stability derivatives were estimated from the equations 
given in the footnote of table I and, are similar to those given in 
reference 5 * 


EESUITS AND DISCUSSION 
Airplane with Wing-lip Tanks 


The results of the calculations to determine the effect of varying 
Cq^ on the neutral-oscillatory— stability boundary of the model with 

wing tip tanks are presented in figure 3 iu the form of a stability chart 
in which the neutral— oscillatory-stability boundaries are plotted as a 
function of the directional-stability derivative C^ and the effective- 

dihedral derivative C 1 . 


The results indicate that, for values of CL 

n P 


of 0 . 25 , 0 . 30 , 


and 0.40, branches of the neutral— oscillatory-stability boundary exist 

for both negative and positive values of C^ , whereas for values of Cq 

P P 

of - 0 . 05 , — 0 . 1 , 0 . 05 , and 0.15 the neutral-oscillatory-stability boundary 
only appears for negative values of C^; that is, in the first quadrant. 

An B = 0 curve for On of — 0.05, -0.1, 0.05 and 0.15 does appear in 


the second quadrant but it is not a neutral— oscillatory-stability boundary 
because the D coefficient for these values of C is negative. Also, 

the dashed portions of the E = 0 curve for C^ of 0.25, 0.30, and 

0.40 in the first quadrant are not neutral-oscillatory— stability boundaries 
because the D coefficient is negative. The relation between the D 
and E boundaries for Cq^ = 0.4 is illustrated in figure 4. The 

coefficient D is negative in the region below the D = 0 boundary. 

For the case where the E = 0 boundary is below the D = 0 boundary, 
there exists a negative D coefficient and hence the dashed portion 
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of the E = 0 curve is- not a neutral— oscillatory— stability "boundary. 
Although, figure 4 is presented only for the case of = 0.40, the 

conclusions drawn from this figure are also applicable to the other 
cases. 


Where two branches of the neutral— oscillatory— stability boundary 
do exist, the significance of the boundaries can be determined by 
investigating the period of the neutr ally stable oscillation for points 
located an each of the branches. As shown in reference 3* the period is 

proportional to i/l therefore , the branch of the neutral— oscillatory- 

stability boundary which is located near the D = 0 boundary, where the 


value of D is small for combinations of 


s 


and C( 




located on 


this branch, represents the neutral— oscillatory— stability boundary of 
a long-period oscillation. 


The part of the curves shown as solid lines for values of 

of 0 . 25 ., 0.30 and 0.40 in the first quadrant (fig. 3 ) represent neutral- 
oscillatory— stability boundaries of a long-period oscillation. The 
other branches of the curves for equal to 0.25, O. 30 , and O.hO 

in the second quadrant^ and the curves of CL equal to -0.05, -0.10, 

0 . 05 , and 0.15 represent neutral— oscillatory— stability boundaries of 
the normal short— period oscillation. The results of figure 3 indicate 

that as C n varies from a negative to a small positive value, a 
P 

stabilizing shift in the oscillatory— stability boundary occurs. Further 
increases in the positive value of C^, however, result in a large 

destabilizing shift of the neutral— oscillatory-stability boundary. 

A better indication of the meaning of the boundaries can be seen 
from a study of the damping and period of the oscillation and the ' 

damping of the aperiodic mode obtained from the roots of the lateral 
stability equation. The reciprocal of the time to damp to one— half 
anqxLitude has been used herein to express damp ing because it was desired 
to express the degree of stability as a direct rather than as an Inverse 
function of the ordinate . Peak positive ordinates consequently indicate 
maximum damping and negative ordinates indicate negative (unstable) 
damping. The results presented in table XX were obtained from roots 
calculated for the case of C^ = 0.40 at 

0.00100 per degree while 


= 0.00227 per degree 
was varied. The conditions 


and CL , _ , 

'dp JT - o - t p 

listed in table H are identified by circled numbers ip figure 3* Table II 
shows that when the neutral— oscillatory— stability boundary is crossed 

(from C, = 0.0025 Per deg to C 


= 0.0025 per deg to 

V 


= P er de S j, the short— period oscil- 


lation becomes stable. This oscillation continues to remain stable as the 
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value of C . 


la is increased negatively. However,, at C 7 =—0.0032 per degree, 
p *3 

a long— period oscillation appears and as the neutral— oscillatory— stability 

"boundary is crossed [from = - 0 . 0 (A -9 per deg to = - 0.00550 per 
degj , the long— period oscillation becoides unstable. Table U also shows 


that the E = 0 curve in the second quadrant 


(s - °- 


is a 


= o.ooiooj 

neutral-oscillatory-stability boundary for the short-period oscillation 
but that the dashed portion of E = 0 curve in the first quadrant is not 
a neutral— oscillatory— stability boundary. 


on 


Besults of calculations made to show the effect of varying C^ 

the period and damping characteristics of the oscillatory find aperiodic 
inodes are shown in figure 5« These calculations were made for a given 
value of CL and CL (point (A) in fig. 3 ) and Bhow that increasing 

p P 

in a positive direction causes an increase in the damping of the 
short— period oscillation throughout the C^ range investigated. 


Airplane without Wing-Tip Tanks 


The effect of C n ^ on the neutral— oscillatory— stability boundary 

for the model without wing-tip tanks is shown in figure 6 . The 
E = 0 curves for values of — 0.05 and —0.10 are neutral— oscillatory- 

stability boundaries for the short-period oscillation. The solid portions 
of the E = 0 curves in the first quadrant for C of 0.10, 0.20. 

n p * 

and 0.40 are neutral— oscillatory-stability boundaries of the lo ng -period 
oscillation, whereas the corresponding branches of the boundaries in 
the second quadrant are neutral-oscillatory-stability boundaries for the 
short— period oscillation. The damping and period data were calculated 
from the roots of the stability equation by varying C n for a given 

value of and C^ (point (B) in fig. 6 ). The results of these 


fc P 


calculations are presented in figure 7 and indicate that the damping cf 
the short-period oscillation increased as CL was increased in a 
positive direction. J? 


A comparison of figures 5 and 7 indicates that shifting fuel from 
the fuselage to wing-tip tankB caused a decrease in the damping of the 
short-period oscillation for all values of investigated. (The 

data of figures 5 and 7 were computed for different c omb inations of 

an( i hut it is assumed that at least a qualitative comparison can 

be made of these results.) 


P 
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For extreme positive Talues of C y,^ in 'both loading conditions the 

damping of the aperiodic modes , "which are normally highly damped, was 
reduced to nearly zero. 


CQHCiosiarrs 


From the results of the investigation made to determine the effect 
an the oscillatory stability of varying (rate of change of yawing- 

moment coefficient with the rolling-angular-velocity factor) the following 
conclusions are drawn: 


1. Increasing 



in a positive direction increased the damping 


of the short— period oscillation for all values of 



2. For any given value of shifting fuel from the fuselage to 

external wing-tip tanks decreased the damping of the short— period 
oscillation. 


Langley Aeronautical Laboratory 

national Advisory Committee for Aeronautics 
Langley Field, Ya., July 23, 19^8 
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CTAR ACTERIBTTC3 OP AIRELAHE USED HI EHB CAICOLAIIOHS 


Items 

Without ving-tip tanks 

With ving-tip tanks 

W, lb 

16,500 

15,700 

w/s 

47-3 

W.9 

b, ft 

39-7 

39-1 

P 

0.00238 

0.00238 

l*b 

15.56 

lk.1 

% 

0.117 

0.302 • 


0.226 

0.1*05 


0.61*3 

0.61*3 

z/b 

0.162 

0.162 

°L 

O.70 

0.70 

a, deg 

12.0 

12.0 

€, deg 

5 

5 

T), deg 

7 

7 

7, leg 

-10 

-10 

P«r radian 

0 . 1*0 

0.1*0 

Oy , per radian 

0 

0 

Cv , per radian 

-0.115 + Cv o/ 

-0.115 + CV 

*P 

^P(tail) 

X P(tail) 

Cjjp, per radian 

- 0 ' 012 + 

-0.012 + CL 

P(tail) 

a 

per radian 

-O.30 + C 7 

T’Ctall) 

-0.35 *t* ^7 

^(tall) 

a C^, per radian 

^Variable + Cn 

^>(tail) 

^Variable + CL 

^P(tail) 

a C 2 ^, par radian 

0.175 + c 7 

^(tail) 

0.175 + Ct_ 

^(tail) 

Cn^,. per radian 

-0.006 + C_ 

■ r (tail) 

-0.006 + Cn 

^(tail) 

CV . per radian 

^P(tail) 

°Variable 

c Variable 


a Tail contributions axe determined from the following equations: 


^(tail) ^ ^(tail) 

C, = 2^- - - sin a) CL 

P(tail) V * ' ^P(tail) 

Cn = c 7 = - 2 sin a)cv 

■^(tail) ^(tail) 13 ^ 15 / (3 (tail) 


^r(tail) 


Kh) °^P(tail) 



^Varied syst emat ic ally as ind ep endent variable to provide the desired range 
of Cj^ for the determination of the stability- boundaries. 

^Varied u _y st emat i c ally as d~nd e p end mr h variable to determine the effect 
of on the stability boundaries. 
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TABLE H 

PERIOD AND RECIPROCAL OF TIME FOR AIRPLANE WITH WING-TIP 
TAMS TO DAMP TO ONE-HALF AMPLITUDE 


Oscillatory mode 


1 


Condition 



''la r T . 

/ , \ (sec) -4/^ 

(per deg) (i/sec) 


^ = 0.00227 



0.0025 

2.10 

-O.O85 

.0010 

2.11 

.183 

-.0032 

2.02 

75-22 

1.048 

.315 

-.00447 

1.96 

17.90 

1.273 

.089 

-.0049 

1.93 

15.99 

1.359 

.003 

- -.0055 

1.90 

14.82 

1.451 

-.089 


CL = 0.0010 


Aperiodic mode 


T. 


1/2 

(1/b6c) 


-0.253 

3.135 

-.216 

2.563 



0.0010 

3.11 

-0.146 

-0.153 

2.793 

-.0015 

3.19 

.861 

-.265 

.900 

-.0030 

2.81 

1.385 

-.575 

.167 

-.oo4o 

2.60 

1.595 

-.870 

.045 
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Figure 1.- System of axes and angular relationship in flight. Arrows indicate 
positive direction of angles. “H = 0 - y - 6 . 







Cfl C»4 






Figure 3.- Effect of the rotary derivative on the neutral -lateral ^oscillatory -stability boundary. 

C'l = 0.70. Airplane with wing -tip tanks. Circled numbers in figure are identified as conditions 
given in table II, The dashed-portion of the R = 0 boundary is not a neutral -lateral-oscillatory- 
stability boundary because the D coefficient of the lateral stability equation is negative. 



D/rechono/-sfc)bil/fy paramete^C 



Figure 4.-- Relative location of the D = 0 boundary wi 
of the R = 0 boundary above the D = 0 boundary r 
Tvyirnrin-rv = 0 40* C'.-r b 0.70. 

— ti - ~ 7 — It “ ' ’ ' ’ 



/ w* r } 


- r\ 
-i/s 


respect to the R = 0 boundary. The portion 
resents the neutral -lateral-oscillatory -stability 
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, Yawing- moment coefficient ctue to 
rotting - ongutar- veto c/fy fa c for , 


Figure 5.- 


Effect of the rotary derivative C n 


on the damping and period of 


the oscillatory and aperiodic modes of a sweptback fighter-type airplane 

with wing -tip tanks. C n = 0.00186 per degree; C 7 = -0.0026 per degree; 

P fr P 


C L =0.70. 



Df reef zonal- sfahif/fy parameter 



Figure 6.- Effect of the rotary derivative on the neutral -lateral -oscillatory -stability boundary. 


A \ n n q txrl tlv\n+ ttH nrr -Hri ton IrCJ 

. u i ^ . j.. . n xuuwuv "" ' ' 


O n(\ HPHo HoaVuaH rvf fViA "R 0 Iv MirviQ.r ^r i o nAt q 

V» I V« J. UVtMUVVt F VJ- ‘’AVili V* A V w WV » '|*»|- ■ * * j AW ■* *W V 


neutral-lateral -oscillatory-stability boundary because the D coefficient of the lateral stability 
equation is negative. 


H 

vo 
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Yawing- /770m on / coofY/c/on/ c/u o t~o 
ro///r?g - argot. lor - ve to c/fy fac tor , 


Figure 7. - Effect of the rotary derivative on the damping and period of 

the oscillatory and aperiodic modes of a sweptback fighter-type airplane 
without wing -tip tanks. C n = 0.0017 per degree; Cj ='0.0017 per 

0 0 


degree; Cl = 0.70. 




